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NITRIC OXIDE IN THE MESOSPHERE AND ITS VARIATIONS 


A. P. Mitra* 

Laboratory for Space Sciences 
NASA Goddard Space Flight Center 
Greenbelt , Maryland 

The distributions, seasonal variations and the resulting 
ionospheric implications of mesospheric nitric oxide are dis- 
cussed. We first consider the principal features of an 
average distribution for this constituent that would be con- 
sistent with photochemical considerations and with the diurnal 
and solar cycle variations in D region ionization and the re- 
lative abundances of the molecular ions NO^ and 0^ . 

We then consider the extent to which the average profile 
may vary seasonally in middle latitudes. Large winter enhance- 
ments occur in NO concentrations at all D region heights above 
about 70 km due to the winter warming at these heights. The 
enhanced region is limited roughly between 70-90 km, the lower 
limit being controlled by the rapid loss of nitric oxide 
through ozone. NO production through the ionic reaction 
O 2 + N 2 -* NO'*’ + NO does not appear to be significant at any 
height. 

Ionospheric implications of these NO distributions are 
also discussed. 


♦NRC-NASA Resident Research Associate, on leave from the National 
Physical Laboratory, New Delhi 12, India. 
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INTRODUCTION 

An important question concerning mesospheric ionization 
and meteorology is the role played by such minor constituents 
as nitric oxide, ozone and atomic oxygen, the metallic species, 
water vapour and ion clusters, and certain excited states of 
the principal constituents. The importance of these in the 
ionization of the mesosphere is now, in most cases, clearly 
established. 

In this work we deal with the distributions, variations 
and the resulting ionospheric implications of only one of these 
minor constituents, eg. nitric oxide. 

Information on the distribution of nitric oxide in the 
mesosphere and the lower thermosphere has been obtained in 
three ways: from rocket measurements of nitric oxide airglow 

in a series of flights by Barth (1966) , from photochemical con- 
siderations (Nicolet, 1965; Barth, 1961; Wagner, 1966), from 
changes in the D region electron density profiles with solar 
activity (Mitra, 1966) and with the time of the day (Mitra, 1968) 
and from pre-sunrise measurements of E region electron density 
(Smith, 1966). In an earlier work (Mitra, 1968) in which estimates 
from these different sources have been reviewed, the following 
empirical distribution was suggested: 

n(NO) = 4xl0~^n (02)e”^^^^'^'^+5xl0 '^n(O) (1) 
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In this work we will obtain some additional ionospheric 
estimates for mesospheric nitric oxide. We find that these 
new values are generally consistent with our previous empirical 
model, escept in the region of minimum aroundSO-lOO km, where 
the new values are larger. We then examine the changes that 
are likely to occur in this distribution in winter from known 
changes in the distribution of atmospheric temperature and in 
the concentrations of atomic oxygen and ozone. The resulting 
concentrations in winter are found to be considerably larger 
than in summer, at heights between 70-90 km, the exact degree 
of enhancement being dependent on latitude. The consequences 
of such changes in the ionization of the winter mesosphere are 
found to be quite severe. 

2. AN AVERAGE DISTRIBUTION FOR NITRIC OXIDE IN THE MESOSPHERE 
AnI) T'hE~i:0WEr THM'O’gPHERE 

Currently available estimates of the nitric oxide con- 
centration in the mesosphere and the lower thermosphere, and 
the sources from which these are derived, are given in Appendix 

1 . 

The unusually large value obtained in Barth's measurement 
has been frequently discussed. In contrast, the different 
ionospheric estimates listed in Appendix 1 are found to be 
generally consistent and in agreement with photochemically 
deduced values . 
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The ionospheric estimates are, however, still far from 
definitive, although these are probably reliable as order-of- 
magnltude values. If the major source in a non-polar latitude 
changes from cosmic rays to La to X rays as one moves from the 
lower to higher mesosphere, then the problem of obtaining a 
reliable value for NO reduces to one of identifying the height 
at which ionization by La predominates. This, it appears, can 
at present be best done from an examination of the diurnal 
variation of electron density profiles. In our previous study 
we used the diurnal profiles obtained by Barrington et al (1963) 
for the high latitude station Kjeller for the period March 17 
to April 7, 1960. This identified the La-dominated levels to 
be around 72 km. A more recent series of results, which we examine 
here in detail, comes from Thrane et al (1968) who obtained de- 
tailed profiles of mesospheric ionization for a large range of 
X values, for both afternoon and morning conditions, from 
partial reflection and cross-modulation experiments conducted 
in the middle latitudes during the IQSY. The afternoon values 

were nearly the same as the morning values (for the same value 

2 

of x) » thus making it possible to use N directly as a measure 

2 

of q, the electron production rate. Plots of N against sec X 
for the heights 65, 70 and 75 km are shown in Fig. 1, along 
with the theoretical variation expected from ionization controlled 
only by La radiation (assuming no diurnal variation for nitric 
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oxide) . We again observe a variation essentially identical with 
that of La at 70 km, a slower variation at 65 km (as would be 
expected if cosmic ray contribution is also important) and a 
faster variation at 75 km, identical with that expected from 
ionization by X rays at 4A° . 

From the results of Fig. 1, we obtain: 

n(N0) s 5x10® cm"^ at 70 km 
8 -3 

n(N0) = 4x10 cm at 65 km (probable value 

-not reliable) 

In all these estimates, as well as in our previous ones 
(Mitra, 1966; Mitra, 1968), the principal concern was to avoid 
any assumption of the loss rate y, given by Y = (1+X) (a^+Xa j^) 
where \ is the ratio of negative ions to electrons, a^^ is the 
ion-ion recombination coefficient and ajj is the coefficient for 
dissociative recombination. This was done because of the current 
uncertainty in our knowledge of Y. Nevertheless, a lower limit 
for the electron production rate can be obtained by writing 
X=0 in the equation: 

" (1+X) (aj)+Xa^)N^ (4) 

If we now limit ourselves to 70 km, for which q is 
negligible, then the minimum value for nitric oxide is set 
by the equation: 
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ai)(NO+)N^ q 
n(NO) = — - ytTT 

where Y(x) = A(NO) Q (La)e defines the ionization 

a 

yield due to La for one NO molecule. 

2 

Estimates of the quantity N /Y(x) from observed electron 

densities at 70 km, over a wide range of solar zenith angles, 

12 —6 

give values in the range (l-3)xl0 cm~ sec, although abnormally 

12 ""6 “1 

large values around 50x10 cm" sec" are also occasionally 

+ 2 

encountered. If q <ar^(N0 )N , and one takes 

C JJ 

ajj(NO^) = 5xl0"^ (T/300)"^ ‘^cm^/s (see also sec. 3), then 

+ —7 3 

ttjjCNO ) at 70 km is around 7x10 cm /s , and 

n(N0) s (0 . 7-2)xl0^cm"^ 

With Mitra's values (1968) of y(= 2xl0"®cm^/s) at 70 km, 

n(N0) w (1 . 4-4)xl0^cm"^ 

A method of considerable value is the use of observations 
of NO^ concentrations, which are now being available with in- 
creasing reliability, along with the laboratory measurements 
of rate coefficients relevant in NO^ chemistry (Mitra, 1968; 
see also Appendix 1) . The reactions relevant in the mesosphere 
are given in Table 1, along with the laboratory values of the 
rate coefficients, where available. It can be shown that NO 
concentration is related to the NO^ concentration in the 
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TABLE 1 

PRINCIPAL REACTIONS FOR NO''" CHEMISTRY IN THE MESOSPHERE AND 

THE LOWER THERMOSPHERE 


NO + hv ^ NO"*^ +e 

+ ^7 + 

Og + Ng 1 > NO + NO 

+ ^5 + 

Oj + NO NO^ + Og 

+ ^11 + 

Ng + 0 NO + N 

"2 + »2 ^ ®2 + "2 

+ ''^1 + 

0^ + Ng i-) NO^ + N 


q(N0) 

i«“17+2» 3 , 

<10 (~10 — )cm /s 

(8+2)xl0"^°cm^/s 

(2 . 5+0 . 8 ) xlO" ^®cm^/s 

(2+l)xl0“^®cm^/s 

1x10 — cm /s 


mesosphere through the relation: 

A- 4 - 

n(N0) [Y (x)+Y 5 n (0*) (0 J)n (N^) (N^) [■ ^^^n{a)^V 9 n( 03 ) <«> 

= ajj(NO'^)n(NO''')Ng 

where the quantity on the right hand side is directly ob- 
tained from observations . 

With the rate coefficients given in Table 3, Eq. (6) takes 
the form: 

A (x ,0g) n (NO) -ttpCNO''') n (N0‘'')Ng- 10"^®n (Ng) n (0+) -B (Og/0) q (Ng) (7) 

with 

A(x,0j) = Y(x)+8xl0"^°n(0j) 
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BCOg/O) = Cl+n(02)/n(0)]"^ 

Since Og is not observed below 80 km, and only exists in negligible 
quantity in the upper mesosphere, the second term on the right 
has no controlling effect in the expression and the uncertainty 
in the value of Vy little effect on the determination NO 
concentration. Direct photoionization of NO ceases to be im- 
portant in NO''^ chemistry above 80 km, but the neutral nitric 
oxide can still control NO^ through process (Og+NO) . Excepting 
for conditions of enhanced nitric oxide, the left hand side is 
negligible above about 100-120 km, and consequently the method 
ceases to be valid above this height. 

Nitric oxide distribution has been computed on the basis 
of Eq. (6) , using NO'*’ distribution obtained by Narcisi and 
Bailey (1965) and with two sets of values for aj-j(N0"^) : 

(a) ajj(N0"’’)=5xl0”^(T/300)”^ *^cm^/s and (b) a constant value 

— 7 3 

of 2x10 cm /s . The distributions so derived are shown in 

Fig. 2 along with other estimates made in this and in previous 

works. For a value of 6 (following Hale, 1968) is used 

for 70 km and 1 for 75 km. The distributions are nearly constant 

6 —3 6 —3 

with concentrations around 6x10 cm~ for (a) and 1.5x10 cm for 

(b) between 70-80 km; the concentrations increase to about 
3xl0^cm ^ at 90 km for (a) and to 7xl0^cm ^ for (b) . 

There is considerable agreement between these values and 
the simplest photochemical calculations. The major reactions 
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TABLE 2 

PRINCIPAL REACTIONS FOR NITRIC OXIDE IN THE MESOSPHERE AND LOWER 

THERMOSPHERE 


MAJOR REACTIONS 


3600+200 


1 . N+0, 


NO+0 


b^=(l . 4^0 .2)xl0”^^e 


3/ 

cm /s 


2 . N+O+M 


la 


> NO+M 


3. N+0 


'lb 


^ NO+hv 


N+NO — ^ N 2 +O 
N+OgC^A) NO+0 


b^^=l.lxl0”^2(^) cm®/s 

bj^^=2.10 ^^cm^/s 
b»=2 .2xl0”^^cm^/s 

b 

“13 “15 3 

k unknown (10 -10~ cm /s used here) 


ADDITIONAL PROCESSES LIKELY TO BE IMPORTANT 


Y 

6. 0J+N2 I.) NO'^+NO 


7 . NO+O3 - 2 -) NO2+O2 

+ ^5 + 

8 . Og+NO — NO +O2 

NO - DISTRIBUTIONS 


Y„=lxlO”^'^i^cm^/s 

_ 1300+100 

b.=9.5xl0"^^e 
4 

Y5=(8_+2)xl0”^^cm^/s 


” 3/ 
cm /s 


1. Mitra’s Empirical Distribution: 

n(N0)=4xl0”^exp(- ^^^- )n(02)+5xl0~^n (0) 

2. Photochemical Equilibrium with 1-4: 

n*(N0) = (0.6+2)xl0'^exp(- 3 . 600-^00 )n(o^) 

+ [(l+0.5)xl0"®+5xl0“^^n(M)^~]n(0) 

3. Photochemical Equilibrium - all processes 1-8: 

n(NQ)-n*(N0) [ ^ n(0)n(N) ] 


1 + 


b^ nW 
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involved in NO photochemistry are given in Table 2 along with 
distributions expected under photochemical equilibrium when the 
major processes are limited to 1-4 and in the more complex case 
including the effect of ionic reactions, possible production 
of NO from the metastable 02(^A) (Hunten and McElroy , 1968), aud 
the very rapid loss of NO through ozone in the lower mesosphere. 
In Figs. 3 and 4 the rates of production of NO (per 
nitrogen atom) and of loss (per NO molecule) are plotted as 
a function of height for heights between 50-120 km. For 
process (Og+Ng) the rate has been normalized by dividing the 
total rate by n(N). For O 2 , N 2 and 0 the distributions used 
are those given by CIRA 1965; for 0^, that given by Maeda and 
Aikin (1967) and for 02(^A), that obtained from the dayglow 
measurements of Evans et al (1968). For atomic nitrogen, three 
distributions are used: N(2) given by Barth (1961) represents 

the lower limit; N(l) obtained by extrapolating downwards the 
photochemical distribution recently derived by Ghosh (1968) 
for heights above 100 km; N(3) representing the upper limit and 
given by Khovistok (1962) (Fig. 5). The temperature-dependent 
NO production process (N+O 2 ) and production through Og(^A) 
compete in importance at heights below 80 km, for k values 
between 10 -lO” cm /s; the uncertainty in the value of k 
makes it difficult to judge which is more important. The 
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former may exceed latter at times of greatly increased temperature 

and for low k-values. Above 80 km, production through process 

(N+0) begins to contribute, with a seasonal variation following 

that of n(0) . The ionic process (Og+Ng) is unlikely to be 

important at any height excepting for very low concentrations 

8 —3 

of atomic nitrogen. If n(N) is 10 cm or larger at these 
heights, as seems more likely, this process is never important. 
Amongst the loss processes (Fig . 4) , the most important is 
(N+NO) for all heights above 70 km (unless atomic nitrogen 
concentration is as low as in Barth's model), but below 70 km 
the process (NO+Og) becomes very effective. 

There is then a sharp cutoff in NO concentration, the cut- 
off level, which depends upon the ratio n(N)/n(0g), is around 
60 km for high atomic nitrogen concentration (model 3) but 
would be as high as 70 km for models 1 and 3. 

6 7 —3 

Since ionospheric estimates yield values around 10 -10 cm” 

8 —3 

at 70 km and (less reliably) 5x10 cm at 65 km, the lower cut- 
off is not acceptable, and the true distribution must at least 
be equal to, and possibly larger than, that of model 1. The 
distribution derived from the complete photochemical equation 

(with the atomic nitrogen distribution of Khovistok and with 
— 14 3 

k=10 cm /s) is shown in Fig. 2 along with the ionospheric 
estimates made in this and the previous works . The agreement 
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Id generally remarkably good except in the region of the minimum 

1 

where the ionospheric estimates are in agreement with the 
photochemical values only when process 5 is included. It has 
also recently been pointed out by Hesstvedt (1968) that when 
one considers the effect of eddy diffusion there is a tendency 
for the minimum to be filled up. Our new estimates confirm this. 

The net result is a nearly constant NO concentration (in the 
neighborhood of 10 cm” ) all the way from 70 to 100 km. 

3. SEASONAL VARIATION IN NITRIC OXIDE 

The very high temperature-dependence of NO production 
through (N+Og) means that any variations in temperature must 
profoundly affect the nitric oxide concentration at levels where 
this process is predominant or competitive, and possibly also 
at other levels through dynamical forces. It is known that 
there is a reversal in the seasonal temperature variation (i.e. 
higher temperature in winter than in summer) from 70 to 90 km, 
at precisely those heights at which this process can become dominant 
with increased temperature. CIRA (1965) gives detailed tables 
of the seasonal variations in temperature at intervals of 10^ 
in latitude. There will thus be large increases in winter 
time nitric oxide concentration at medium and high latitudes. 

In a static atmosphere, this, in turn, would produce large 
localized regions of enhanced electron density at levels where 
the (N+Og) process contributes. 
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The summer-to-winter changes in nitric oxide as a result 

of the increase in mesospheric temperature are given in Fig. 6 

for two representative latitudes 50^N and 60^N, and represent 

conditions on January 1 and July 1, NO production through metastable 

02 ^ A) is included and its coefficient is assumed to have a value 
-14 3 

of 10 cm /s, independent of temperature. If the loss of NO through 

process (NO+O^) is ignored, one finds a decrease in nitric 

8 —3 

oxide concentration from summer to winter from 6x10 cm to 

SxlO^cm ^ at 60 km at 60^N latitude and from 3.5xl0^cm ^ to 
7 —3 o 

5.5x10 cm at 50 N, the decrease being roughly a factor of 
10 at all heights below 60 km. Above about 67 km, however, 
the situation is drastically different. At these heights there 
is a rapid increase in NO concentration; the increase is almost 
1 order of magnitude at 75 km, and similar at 80 km. This large 
increase in NO concentration at altitudes above 67 km continues 
up to about 90 km. The entire seasonal variation in NO con- 
centration for 50°N is shown in Fig. 7, The large enhancement 
in NO concentration from October to March has a striking resemblance 
to HF absorption for similar latitude areas. 

4. IONOSPHERIC IMPLICATIONS OF WINTER ENHANCEMENT IN NITRIC 
OXIDE 

These large enhancements in NO concentration at medium and 
high latitudes has serious implications in D region ionization 

and especially in its seasonal variation. At 75 km and at 

o -2 -3 

50 N, q(N0) jumps from a value of 4x10 cm sec 


in July to 
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—2 —3 —1 1 

about 30x10 cm sec in January if production through 02( A) 

-2 -3 -1 

is ignored, and from 30x10 cm sec in July to about 

_2 -3 -1 ] 

38x10 cm sec in January if 02< A) process (with a rate 

— 14 3 — L 

coefficient of 10 cm sec ; is included. As will be seen 

in Fig. 7, q(N0) is large over the entire period October to 

March. At 60°N, with 02(^A) process, q(N0) has a nearly constant 
-2 -3 -1 

value of 20x10 cm sec from February to November, with some 
increase at the equinoxes and somewhat lower values in January 
and December . 

In the summer mesosphere and when NO production through 

02(^A) is negligible, the X ray contribution exceeds the 

contribution due to La at 80 km, during all phases of solar 

activity, if NO concentration is not much in excess of that 

given by curve 1 in Fig. 2. With a predominating X-ray source, 

the winter production rate would then be reduced by a factor of 

about 5x10 ^ at 50^N. However, the large enhancement in nitric 

oxide concentration in winter, yielding a production rate of 
-1 -3 -1 

8x10 cm sec in January entirely reverses the situation. The 

total rate of electron production at 80 km, as a result of this, 

— 1 —3 —1 —2 —3 —1 

becomes 8x10 cm sec in winter against 6x10 cm sec in 

summer. One would thus expect a four fold increase in electron 

density in winter from the enhancement of nitric oxide alone. 

When 02( A) process is included (with k=10 cm /s and assumed 

independent of temperature). La is found to dominate both summer 

and winter ionization at 80 km, but winter production rate 

continues to be high, although now only twice as high as in summer. 

At lower altitudes, the wintertime enhancement in nitric oxide 
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becomes progressively smaller and consequently the winter 
electron density will approach the summer value and will 
ultimately fall below it. 


TABLE 3 

WINTER ANOMALY AT 50°N DUE TO NITRIC OXIDE ENHANCEMENT 



SUMMER (27*^) 

WINTER (73°) 

q(La) 

8xl0~^(4.6xl0"^^^ 

“X *“3 ™ X 

8x10 cm” sec” (10xI0~^^ 

q(X) 

5x10“^ 

o in-4 -3 -1 

2x10 cm sec 

q(total) 

6x10”^(5x10“^^^) 

8x 10” ^cm”^s ec” Vox ID”^^^) 

N ♦ 
e 

3 . 5xl0^(lxl0^^^) 

3xl0^cm”^(l .4xl0^^^) 


— 7 3 

*Y=5x10 cm /s ♦♦Including 0 (^A) process 

2 


Two recent experimental observations, relevant to the 
above discussion, are (i) the electron density profiles of 
Gregory (1968) obtained with partial reflection at 44°S location 
and (ii) the observations made by Ferraro and Lee (1967) at the 
Pennsylvania State University (40°N) during the IQSY with the 
cross modulation experiment. The Penn State results, although 
showing considerable scatter, show no evidence of winter anomaly, 
and is consistent with the variation expected from ionization 
by X-rays alone. Gregory's observations showing winter electron 
densities as large as in summer would be consistent with a 
wintertime enhancement of the type mentioned above, but by 
a smaller magnitude. 
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It is also important to note that the winter enhancement 
in nitric oxide is limited to heights above 70 km, below which 
the winter mesosphere is cooler than the summer mesosphere. 

The sharp decrease in NO concentration due to the rapid loss of 
NO through Og can raise this level to about 75 km in case of 
low atomic nitrogen concentration. Present evidence indicating 
that the winter anomaly is essentially confined to height between 
80-90 km (Gregory, 1968; Sechrist et al, 1968) is in agreement 
with this . 
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Fig. 1 


Fig. 2 


Fig. 3 


Fig. 


LEGEND OF FIGURES 

2 

- Plots of N (proportional to the electron production 
rate) against log sec X at selected heights of 65 , 

70 and 75 km. Electron density observations are 

those by Thrane et al (1968) . Solid lines and theoretical 
variations of q(La) with sec X, and the dotted line 
for 75 km is the variation of the electron production 
rate due to 4A alone. 

- Different ionospheric estimates of mesospheric nitric 
oxide, along with two photochemical distributions 
one with and the other without the 02 (^A) process. 

The dotted portion of the theoretical curves shows 
the effect of rapid loss of NO through 0^ for the 
atomic nitrogen distribution of Khostikov. 

- Relative magnitudes of the specific rate of production 

of NO through different processes (expressed per 

nitrogen atom). For the ionic reaction O 2 +N 2 -* NO^+NO , 

— 15 3 

an upper limit of 10 cm /s has been used for the rate 
coefficient. The two curves refer to the atomic nitrogen 
models N(l) and N(2) in Fig. 4. The corresponding 
curve for N(3) , will be smaller than that for N(l). 

- Assumed distributions of atomic nitrogen and ozone in the 
mesosphere . 
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Fig. 5 


Fig . 6 


Fig. 7 


- Relative magnitudes of the NO loss rates for different 
reaction, (expressed per NO molecule). The two curves 
for (N+NO) refer to the atomic nitrogen distributions 
N(l) and N(2) in Fig. 4. 

- Calculated NO distributions for January and July in 
the mesosphere corresponding to CIRA temperatures. 

The sharp cut-offs at the lower levels, different 
for different atomic nitrogen distributions, should 
be noted. Production of NO through [N+OgC^A)] is 
included with constant rate coefficient of 10”^”^cm^/s. 

- Seasonal variation of nitric oxide and the resulting 

electron production rates calculated for 50°N and for 

a level of 75 km. Calculation are given for: (a) 

N+ 02 (^A) process negligible and (b) a temperature- 

— 14 3 

independent rate of 10 cm /s for the process 
N+OgC^A) - NO+0. 
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